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ABSTRACT 

Using mean relative peculiar velocity measurements for pairs of galaxies, we 
estimate the cosmological density parameter fl m and the amplitude of density 
fluctuations <r 8 . Our results suggest that our statistic is a robust and reproducible 
measure of the mean pairwise velocity and thereby the Q m parameter. We get 
VL m = 0.30^o;q7 and cr 8 = LIS^; 2 . 2 . These estimates do not depend on prior 
assumptions on the adiabaticity of the initial density fluctuations, the ionization 
history, or the values of other cosmological parameters. 

Subject headings: peculiar velocities, the cosmological density parameter 

1. Introduction 

In this paper we report the culmination of a program to study cosmic flows. In series 
of recent papers we introduced a new dynamical estimator of the Q m parameter, the dimen- 

1 Dept. of Physics & Astronomy, University of Kansas, Lawrence, KS 66045, USA 

2 Racah Institute of Physics, Hebrew University, Jerusalem, Israel 

3 Institute of Astronomy, Zielona Gora University, Zielona Gora, Poland 

4 Copernicus Astronomical Center, Warsaw, Poland 

5 Observatoire de la Cote dAzur, Nice, France 

6 Astrophysics, University of Oxford, Oxford 0X1 3RH, UK 

7 Dept. of Astronomy, University of California, Berkeley, CA 94720-3411, USA 

8 IEEC/CSIC, 08034 Barcelona, Spain, and INAOE, Astrofisica, Puebla 7200, Mexico 

9 Dept. of Physics, University of Florida, Gainesville FL 32611-8440, USA 
10 Blackett Laboratory, Imperial College, London SW7 2AZ, UK 
11 European Southern Observatory, D-85748 Garching, Germany 
12 Physics Dept., Carnegie Mellon University, Pittsburgh, PA 15213, USA 
13 Center for Radiophysics and Space Research, Cornell University, Ithaca, NY 14853, USA 
14 Dept. of Physics & Astronomy, Dartmouth College, Hanover, NH 03755-3528, USA 



- 2 - 



sionless density of the nonrelativistic matter in the universe. We use the so called streaming 
velocity, or the mean relative peculiar velocity for galaxy pairs, i>i2(?"), where r is the pair 
separation (Peebles 1980). It is measured directly from peculiar velocity surveys, without the 
noise-generating spatial differentiation, used in reconstruction schemes, like POTENT (see 
Courteau et al. 2000 and references therein). In the first paper of the series (Juszkiewicz et 
al. 1999), we derived an equation, relating i>i2(r) to Q m and the two-point correlation func- 
tion of mass density fluctuations, £,{r). Then, we showed that v\2 and Q m can be estimated 
from mock velocity surveys (Ferreira et al. 1999), and finally, from real data: the Mark III 
survey (Juszkiewicz et al. 2000). Whenever a new statistic is introduced, it is of particular 
importance that it passes the test of reproducibility. Our Mark III results pass these tests: 
the v\2(r) measurements are galaxy morphology- and distance indicator-independent. 

In this Letter we extend our analysis to three new surveys, with the aim of testing 
reproducibility on a larger sample and, in case of a positive outcome, improving on the 
accuracy of our earlier measurements of Q m and <t 8 , the root-mean-square mass density 
contrast in a sphere of radius of 8 h^ 1 Mpc, where h is the usual Hubble parameter, Hq , 
expressed in units of 100 kms _1 Mpc _1 . In our notation, the symbol as always refers to 
matter density, while <7g SCz refers to the number- density of PSCz galaxies. 

Unlike our analysis, other estimators of cosmological parameters are often degenerate, 
hence o~ 8 and Q m can not be extracted without making additional Bayesian prior assumptions, 
which we call conventional priors: a particular choice of values for h, the baryon and vacuum 
densities, and Q\, the character of the primeval inhomogeneities (adiabaticity, spectral 
slope, t/s ratio), the ionization history, etc. (Bridle et al. 2003). The estimates of Q m and 
(Tg presented here do not depend on conventional priors. The only prior assumption we make 
is that up to (T 8 , the PSCz estimate of £(r) describes the mass correlation function. We test 
this assumption by comparing the predicted Vu(r) to direct observations. We also check 
how robust our approach is by replacing the PSCz estimate of £(r) with an APM estimate 
and two other pure power-law toy models. 



2. The pairwise motions and galaxy clustering 

The approximate solution of the pair conservation equation derived by Juszkiewicz et 
al. (1999) is given by 

Wl2 (r) = -|iJ r^ 6 |(r)[l + a|(r)], (1) 
= , s 3 f,T £ (x) x 2 dx 

where a = 1.2 — 0.657, an d 7 = — (dln^/dlnr)|^=i. As a model for £(r), we use the Fourier 
transform of the PSCz power spectrum (Hamilton & Tegmark 2002, eq.[39]), which can be 
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expressed as 

£(r) = (a 8 /0.83) 2 [ (r/n)^ + (r/r 2 )"^ ] , (3) 

where 77 = 2.33 h~ x Mpc, r 2 = 3.51 h~ x Mpc, 71 = 1.72, 72 = 1.28, and 0% is a free parameter. 
If the PSCz galaxies follow the mass distribution, then <7g = <jg SCz = 0.83. The quantities 
erg and £(r) describe nonlinear matter density fluctuations at redshift zero. The PSCz fit 
with <7 8 = 0.83 in eq. (3) is plotted in figure 1, together with the APM correlation function 
measurements for comparison. For r < 15 h~ x Mpc, the APM correlation function is well 
approximated by eq. [3] with 77 = 3.0 hr x Mpc, r 2 = 2.5 h~ x Mpc, 71 = 1.9 and 72 = 1.1. For 
2 h^ 1 Mpc < r < 15 h^ 1 Mpc, which is the range of separations of interest here, the PSCz 
and APM correlation functions in figure 1 are almost indistinguishable. This provides an 
added reason to believe that choosing PSCz as a template for £(r) was a good idea. To test 
the stability of our conclusions with respect to uncertainties regarding the small-r behavior 
of £(r), we will compare predictions for i>i 2 (r) based on PSCz parameters for eq. [3] with 
those based on the APM survey. To study the sensitivity of i>i 2 (?") and inferred cosmological 
parameters to the assumed slope of £(r), we will also consider two simplified, pure power-law 
toy models, given by 

£(r) = (ag/0.83) 2 (r/r H , (4) 
where 7 = 1.3 and 1.8, while r = 4.76 h~ l Mpc and 4.6 h^ 1 Mpc, respectively. 



3. Peculiar velocity surveys 

We will now describe our measurements. Each redshift-distance survey provides galaxy 
positions, f A , and their radial peculiar velocities, s A = r A -v A /r A = f A -v A , rather than three- 
dimensional velocities v A . We use hats to denote unit vectors while indices A, B = 1,2,... 
count galaxies in the catalogue. Consider a set of pairs (A, B) at fixed separation r = \r AB \ , 
where r AB = f A — r B . To relate the mean radial velocity difference of a given pair to 172(7), 
we have to take into account a trigonometric weighting factor, 

{sa-Sb) = v 12 (r)q AB q AB = f AB -(f A + f B )/2 = -q BA . (5) 

To estimate v 12, we minimize the quantity x 2 { v i2) = ^2 A B [{ s a — s B ) — q AB Vi2(r) ] 2 . The 
condition dx 2 / dv± 2 = implies 

vu(r) = ^2(s A - s B ) q AB / ^2q AB 2 . (6) 

A,B A,B 

In this study we use following independent proper distance catalogues. 

1. Mark III. This survey (Willick et al. 1995, 96, 97) contains five different types of 
data files: Basic Observational and Catalogue Data; Individual Galaxy Tully-Fisher (TF) 
and D n -a Distances; Grouped Spiral Galaxy TF Distances; and Elliptical Galaxy Distances 
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as in the Mark II (for TF and D n -a methods, see Binney & Merrifield 1998). The subset we 
use here contains 2437 spiral galaxies with TF distance estimates. The total survey depth 
is over 120 hr 1 Mpc, with homogeneous sky coverage up to 30 hr x Mpc. 

2. SFI (da Costa et al. 1996; Giovanelli et al. 1998; Haynes et al. 1999a,b). This is an 
all-sky survey, containing 1300 late type spiral galaxies with /-band TF distance estimates. 
The SFI catalogue, though sparser than Mark III in certain places, covers more uniformly 
the volume out to 70 hr 1 Mpc. 

3. ENEAR (da Costa et al. 2000). This sample contains 1359 early type elliptical 
galaxies brighter than m B = 14.5 with D n -a measured distances. ENEAR is a uniform, 
all-sky survey, probing a volume comparable to the SFI survey. 

4. RFGC (Karachentsev et al. 2000). This catalogue provides a list of radial velocities, 
HI line widths, TF distances and peculiar velocities of 1327 spiral galaxies that was compiled 
from observations of flat galaxies from FGC (Karachentsev et al. 1993) performed with the 
305 m telescope at Arecibo (Giovanelli et al. 1997). The observations are confined within 
the zone 0° < 5 < +38° accessible to the radio telescope. 

Figure 2 shows our estimates of Vi2(r). Although the catalogues we used are indepen- 
dent, distinct and survey very different galaxy and morphology types, as well as different 
volumes and geometries, our results are robust and consistent with each other. The error 
bars are the estimated 1-a uncertainties in the measurement due to lognormal distance er- 
rors (around 15%), sparse sampling (shot noise), and finite volume of the sample (cosmic 
variance). For more details on error estimates used here, see Landy & Szalay (1992), Haynes 
et al. (1999a,b), and Ferreira et al. (1999). 

The agreement among the V\2{r) estimates from different surveys, plotted in figure 2 
becomes even more impressive when compared to discrepancies between different estimates of 
a close cousin of our statistic, the pairwise velocity dispersion, <T 12 (r). The velocity dispersion 
appears to be less sensitive to the value of Q m than to the presence of rare, rich clusters in 
the catalogue and to galaxy morphology, with estimates of au at separations from one to a 
few Mpc varying from 300 to 800 km s -1 from one survey to another (Davis & Peebles 1983; 
Zurek et al. 1994; Marzke et al. 1995; Zhao et al. 2002). The lack of systematic differences 
between Vu(r) estimates in figure 2 is incompatible with the linear biasing theory unless 
the relative elliptical-to-spiral bias, &e/&s> is close to unity at separations r > 5/i~ 1 Mpc, in 
agreement with our earlier studies (Juszkiewicz et al. 2000); for the same reason our results 
strongly disagree with recent semi-analytic simulations (Sheth et al. 2001; Yoshikawa et 
al. 2003). 

In figure 3 we show the results for each of the catalogues we investigated, as in figure 2, 
but now we overlay the weighted mean of the individual catalogues. Since the results are 
robust, combining the catalogues reduces the errors and gives us a strong prediction for the 
parameter values. Figure 3 shows the results of our theoretical best fits: the solid (dotted) 
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line follows the double power law correlation function using the PSCz (APM) correlation 
function (eq. [3]). Clearly, the slope differences in £(r) at small separations do not affect 
t>i2(V) in the range of separations we consider. Moreover, given the error bars on Vi 2 , the 
7 = 1.3 power-law toy model prediction for Vi 2 (r), as well as the resulting best fit values 
of a 8 and Vl m are similar to those based on the APM and PSCz correlation functions. For 
a 8 PS 1 and £(r) oc r -7 at r > 10 hr x Mpc, linear theory applies and V\ 2 oc r 1-7 . Therefore all 
three of the models considered above give v\ 2 oc r~°' 3 , in good agreement with the observed 
nearly flat Vi 2 (r) curve. All of the above does not apply to our 7 = 1.8 toy model, which 
is significantly steeper than the APM and PSCz £(r) at large r, and for a 8 1, the Vi 2 (r) 
is expected to drop almost by half between 10 and 20 Yc 1 Mpc. It is possible to flatten the 
v\ 2 {f) curve only by increasing a 8 and extending the nonlinear regime to larger separations. 
The example considered here gives a 8 — 1-76, in conflict with all other estimates of this 
parameter (see the discussion below). Correlation functions, steeper than APM or PSCz 
often appear in semi-analytic simulations and this example shows how Vi 2 (r) measurements 
can be used to constrain those models. 

In figure 4 we plot the resulting 1,2,3 and 4a likelihood contours in the (fl m , 08 ) plane. 
The quoted errors define the la, or 68% statistical significance ranges in each of the two 
parameters and correspond to the innermost contour in figure 4. The low % 2 per degree of 
freedom is indicative of the correlations between Vi 2 (r) measurements at different separations 
r. One of the sources of correlations is the finite depth of our surveys. Note also that since 
we are dealing with pairs of galaxies, the same galaxy can in principle influence all separation 
bins. The contours derived using the PSCz correlation function (eq. [3]), are shown in the 
bottom right panel. The best fit values are 

n m = 0.30l£S? and a 8 = 1.13l°;| . (7) 

The likelihood contours based on the APM correlation function (with best fit values Q m = 
0.34+°;^ and a 8 = 1-151°;^) and the 7 = L3 power-law model (fl m = 0.23+°;^ and 
erg = 1. 20^^25) are similar. Our estimate of a 8 agrees with the results of studies of clustering 
of galaxy triplets in real and Fourier space in three different surveys: the APM (Gaztanaga 
1995; Frieman & Gaztanaga 1999), the PSCz (Feldman et al. 2001) and the 2dF (Verde 
et al. 2002). A similar value of a 8 was recently inferred from the observed position of 
the inflection point in the APM £(r) (Gaztanaga & Juszkiewicz 2001). All of the above 
measurements are consistent with a a 8 within 20% of unity. A a 8 close to unity follows from 
maximum likelihood analysis of weak gravitational lensing (Van Waerbeke et al. 2002) after 
assuming £l\ = 0.7, Vl m = 0.3 and h = 0.7. Measurements of the abundance of clusters 
(Bahcall et al. 2003) tend to give a 8 closer to the lower end of our 68% interval if Q m = 0.3. 
The good agreement between these results, obtained with different methods, riddled with 
systematic errors of different nature, suggests that our estimates of statistical errors are 
reasonable and that the systematic errors are subdominant (unless there is an evil cosmic 
conspiracy of errors). The parameters in eq. [7] also agree with those inferred from the 
power spectrum of the anisotropy of the cosmic microwave background (CMB) temperature 
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distribution on the sky: a 8 = 0.9 ± 0.1, and Vt m = 0.29 ± 0.07 (see Table 2 in Spergel 
et al. 2003). It is important to bear in mind, however, that unlike the CMB results, our 
estimates were obtained from the velocity and PSCz data alone, without the conventional 
priors. Therefore the Vu(r) measurements combined with the CMB or the supernova data 
can be used to break the cosmological parameter degeneracy. Choosing 7 = 1.8, which is 
significantly steeper than the observed £(r), gives Vt m = 0.141q;q4 and erg = 1.76ta26 (hg ure 
4, upper right), in conflict with all of the independent estimates of as, discussed above. This 
suggests that the observed slope of the APM and PSCz correlation functions is close to the 
slope of the dark matter correlation function. 
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Fig. 1. — The APM correlation function measurements (circles with error bars) compared to 
four closed-form expressions for £(r): two power-law toy models with slopes 7 = 1.3 (short- 
dashed line) and 1.8 (long-dashed line) and two more realistic, broken power-law empirical 
fits, given by eq. [3]. The latter two represent the PSCz (solid line) and the APM survey 
(dotted line). All four expressions for £(r) assume as = 0.83. 
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Fig. 2. — The pairwise velocities Vi2(r) for the four surveys. The Mark III-S Vi2(r) measure- 
ments come from our earlier work (Juszkiewicz et al. 2000). Clearly, the results from all 
surveys agree well with each other. 
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Fig. 3. — The crosses and the associated error bars show the weighted mean pairwise velocity, 
obtained by averaging over four surveys. Individual survey data points are also shown; we 
have suppressed their error bars for clarity These direct measurements of vu are compared 
to four vu^r) curves, derived by assuming four different models of £(r), plotted in figure 1. 
The labels identify best fit fl m and <r 8 parameters. 
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Fig. 4. — The results of the maximum likelihood analysis. The upper panels show results 
for power-law toy models, while the bottom panels are based on realistic representations of 
observations: the APM and PSCz data, respectively. Likelihood peak coordinates and the 
values of \ 2 for each model are also indicated. The innermost contours define the 68%, or 1-a 
areas around the peaks. The remaining nested contours show the 2, 3 and 4-cr boundaries. 



